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Compaction dynamics of metallic nano-foams: A molecular dynamics simulation study 
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We investigate, by molecular dynamics simulation, the generic features associated with the dy- 
namic compaction of metallic nano-foams at very high strain rates. A universal feature of the 
dynamic compaction process is revealed as composed of two distinct regions: a growing crushed 
region and a leading fluid precursor. The crushed region has a density lower than the solid material 
and gradually grows thicker in time by snowplowing. The trapped fluid precursor is created by 
ablation and/or melting of the foam filaments and the subsequent confinement of the hot atoms in 
a region comparable to the filament length of the foam. Quantitative characterization of nano-foam 
compaction dynamics is presented and the compacted form equation-of-state is discussed. We argue 
that high-energy foam crushing is not a shock phenomenon even though both share the snowplow 
feature. 

PACS numbers: 52.35.Py, 52.35.Qz 
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INTRODUCTION 

The emergence of high-energy laser systems, such as 
the National Ignition Facility (NIF) , presents the oppor- 
tunity to create and interrogate states of matter under 
unprecedented extreme conditions of pressure, temper- 
ature, and strain rates. NIF experiments will also be 
used to conduct astrophysics and basic science research 
and to develop carbon-free, limitless fusion energy. In the 
NIF Energy program, a metallic nano-foam shell has been 
proposed as an important ingredient in the design of the 
double-shell Inertial Confinement Fusion (ICF) target. 
The dynamics of its collapse upon implosion wiW play 
an essential role in achieving inertial confinement fusion. 
However, what is lacking is a fundamental understanding 
of the relationship between the microphysics of the ma- 
terials properties of the drive impactor ( "reservoir" ) and 
the target materials. An atomic level simulation effort 
can serve as the basis to relate the microstructurc of the 
reservoir impactor to the pressure profile on the target. 
The results of these ultra-scale atomic level simulations 
can be significant in improving the predictive capabil- 
ities of continuum level simulations shock compression 
NIF experiments. 

The study of matter under planetary interior condi- 
tions of high pressure and relatively low temperature en- 
vironments is an important research area for NIF 
In the traditional laser-based EOS experiment, a strong 
shock is launched in a material which instantaneously 
increases both the temperature and pressure tracing the 
Hugoniot line in the phase diagram. However, broader 
knowledge of the equation-of-state is required. Recent 
work by Remington and co-workers [2tJ| has demon- 
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strated the possibility of a quasi-isentropic compression 
on a laser based platform. The intense laser pulse drives 
a reservoir which unloads across a gap onto the sample 
of interest. Low density porous materials are used as a 
reservoir because they provide a convenient method to 
independently vary density without changing the atomic 
number. The basic idea is to vaporize the foam, and com- 
press the hot vapor, thereby increasing the pressure and 
temperature of the "vapor piston" smoothly. Molecular 
dynamics simulation techniques can provide design tools 
to tailor the pressure pulse for isentropic compression of 
a solid. 

In this study, our simulation systems are compos- 
ite crystal/nano-foam sandwiches. We have constructed 
physical models and developed the simulation and anal- 
ysis programs with supporting visualization tools for ap- 
plication to study high strain rate compaction of metal- 
lic nano-foam by large-scale computation. Our study 
incorporates computer-generated nano-foams created by 
rapidly quenching a high temperature, phase-separating 
fiuid. Our computer foam has been used earlier to study 
surface-stress-induced relaxation of Au nanostructures 
5] . This present study elucidates the generic features as- 
sociated with the dynamic compaction of metallic nano- 
foams at very high strain rates. Quantitative character- 
ization of nano-foam compaction dynamics is presented 
and equation-of-state of the compacted form is discussed. 



COMPUTER MODELING FOR THE 
SIMULATION 

Our simulation method is Molecular Dynamics (MD), 
where it is assumed that the motion of the atoms is gov- 
ern by Newton's second law 0, 0] . The next ingredient 
is the chosen interatomic potential. It can be simple or 
complicated. While choosing a complicated potential to 
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describe a particular material is often desired, the physics 
of a complicated process can be more transparent and 
discovery of generic behavior may be more readily forth- 
coming with the choice of a simple potential. We assume 
the Voter-Chen EAM potential for copper as our simple 
potential [1] . We acknowledge the limitations associated 
with this choice of potential and will discuss some, such 
as the neglect of ionization. However, we emphasize that 
our goal is to elucidate the generic features of the com- 
paction dynamics of nano-foams. 



nally, the re-scaled distance field was used to carve out 
atoms with a pore filament topology that was identical to 
the original sample, but with specified scales, density pro- 
files and atomic arrangements 13, These steps are 
illustrated in Fig. [5] Further details are to be published 
and a movie of the nano-foam structure may be viewed 
at 'http: / /www.llnl.gov/largevis /atoms /challenge2007| . 
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FIG. 1: Making computer foam by spinodal decomposition. 
The method is the rapid quenching of a superheated fluid into 
the two-phase region followed by the freezing in early time of 
the phase-separating liquid. 

Making computer foam was our first step. Wc took 
the copper atoms to their super critical state and rapidly 
quenched it to the middle of its two-phase fiuid region of 
its phase diagram where it underwent phase separation 
by spinodal decomposition [9] . See Fig. [TJ We stopped 
the fiuid phase separation process by rapidly quenching 
to the solid region, giving rise to a nano-porous foam. 
The initial porous sample was small and the filaments 
were amorphous. 

For studies of the material response of pore structures 
under various conditions, it is desirable to construct pore 
structures at larger length scales, with controlled atomic 
arrangements (e.g. perfect crystal, grains, etc.), and with 
specified density profiles or filament cross-sections. Be- 
low, we briefiy outline a topological filtering approach to 
analyze and synthesize such designed pore structures us- 
ing the original small, amorphous sample as a template. 

First, an initial set of atomic positions from the spin- 
odal decomposition process was recorded. Next, a prox- 
imity field of original atoms on a regular grid was made. 
Then a signed distance field relative to the solid/ void in- 
terface surface was computed. A topology preserving sur- 
face propagation was performed from the interface sur- 
face to produce a topologically clean distance field and 
curved skeleton of the pore filament structure. A dis- 
tance field was computed from the curved skeleton and 
then re-scaled to produce uniform density profiles. Fi- 



FIG. 2: Graphically outlining the topological filtering ap- 
proach for analyzing and synthesizing nano-pore structures. 
See text for explanation. 

For the three-dimensional pore structure used in this 
study, a solid/void interface was determined (shown as a 
green surface in the movie) , along with the solid filament 
center lines (shown in red in the movie). We have created 
foams with filament sizes up to 35 nm while experiments 
are creating foams with filaments of solid copper between 
20 and 100 nm in diameter. 

We used this computer foam to construct the system 
for the shock simulations. We call it a sandwich, and 
its building blocks are comprised of foam, crystal and 
vacuum constructs. We will study the sandwich configu- 
ration shown in Fig. [31 
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FIG. 3: The simulated sandwich configuration. Its building 
blocks are comprised of foam, vacuum and crystal constructs. 

The purpose of the nano-porous foam block is to pro- 
vide a source for creating hot vapor by high-energy im- 
pact. The anticipated history can be easily summarized. 
The piston moves at constant velocity to the right, crush- 
ing and heating the foam block by a rapidly moving front. 
We describe this as a compaction front in contrast to 
the popular shock front description. It was anticipated 
that at sufficiently high piston velocities, the foam com- 
pletely vaporizes. Hot vapor expands from the heated 
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foam into the vacuum and is compressed by the driving 
piston. This, in turn, creates a continually rising pres- 
sure on the face of the right-hand target crystal. In this 
paper, our primary interest is the qualitative features of 
the equation of state of the compacting foam. The initial 
density of the computer foam in this study is taken to be 
15% of the density of solid copper. 

THE SIMULATIONS 
Universal dynamical features of a shocked sandwich 

A universal feature of the compaction process traveling 
through the porous solid has emerged from our simula- 
tions. The feature has two regions: the growing com- 
pacted region and the trapped fluid precursor. See Fig. U) 





distance (arbitrary unit) 

FIG. 5: Growing densification of the crushed foam as a con- 
sequence of snowplowing. The vertical dash lines define the 
piston position and the solid lines are the respective instanta- 
neous density profiles. Superimposed in the background are 
snapshots of the sandwich. 
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FIG. 4: Top: the nano-foam prior to impact. Bottom: the 
universal feature of the compaction process: the compacted 
foam (at left) and the vapor precursor (at right). The colors 
represent the potential energy of the individual atoms. The 
limits are blue being the energy of a solid copper atom and 
red being a weakly bonded or isolated atom. 

In contrast to a shocked solid which achieves con- 
stant density immediate behind a sharp front, foam com- 
paction is an evolving phenomenon where growing densi- 
fication occurs as the region falls back from the forward 
moving broad front as shown in Fig. [5l The compacted 
region has a peak density lower than its cold solid mate- 
rial. A feature shared by shocking and by compacting is 
that growth is achieved by snowplowing. 

A trapped fluid precursor is created by ablation and/or 
melting of the foam filaments and the subsequent trap- 
ping of the hot atoms within a mean-free-path that is 
comparable to the filament length of the foam. The 
high-energy fluid atoms do not experience free expan- 
sion. Instead, they traverse a maze of uncorrelated, ran- 
dom channels restricting free flow to the vacuum. This 
fluid precursor can be liquid or vapor, depending on the 
piston impact speed. See Fig. [B] The fluid precursor 



FIG. 6: The compaction precursor breaking out into the vac- 
uum. This fiuid precursor can be liquid or vapor, depending 
on the piston impact speed. The individual snapshot times 
are normalized so that the apparent rates of collapse of the 
four sandwiches are the same. See caption of FigO for de- 
scription of colors. 



rapidly achieves an approximately fixed length, travel- 
ing at a speed comparable to the compaction front, Uf. 
The compaction front speed is greater than the piston 
speed. Up, but slower than the shock speed in the perfect 
solid. This picture is consistent with a recent experiment 
where a calorimeter measurement suggested that a sig- 
nificant temperature rise preceded the compaction front 
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traveling through a foam [1/ 



Features of the ejecta leaving the foam 

We drove our sandwich with pistons that travelled at 
speeds of interest to recent experiments and discovered 
a phase change associated with the ejecta at the slower 
speeds. This is best seen as the ejecta leaves the foam 
boundary as shown in Figures [6] and [T] 

For Up < 1 km/s, there was only foam compaction and 
no ejecta. As Up was increased from 1 km/s to 4 km/s, we 
observed a liquid front, passing to liquid jetting, liquid 
cluster spraying, and, at the highest speed, hot vapor 
ejecta. Since the recent EOS experiments are presently 
approaching Up ~ 20 km/s, we can conclude that the 
ejecta will be a hot vapor at these higher piston speeds. 



1.0 




FIG. 7: Snapshots of the liquid ejecta at Up = 1.6 km/s. See 
caption of Fig[5]for description of colors. 



EOS simulations of the compacting foam 

Figure [8] shows the respective density profiles for pis- 
ton speeds from 1 to 20 km/s at decreasing times so as 
to capture the same piston location. The chosen piston 
location is where the compacted regions are approach- 
ing their greatest thickness; i.e., the respective fronts are 
near breakout or the foam- vacuum interface. 

Figure [5] also shows that the density profiles are nearly 
identical for Up > 5 km/s. There is not perfect overlap of 
the density profiles in the region around 3, 500 A, but the 
variation is not correlated with Up. The commonality of 
the different profiles is striking and quite unexpected. 

Extending the Up ^ 5 km/s simulation to four times 
in time and distance, we note that the smoothed profile 
of the compacted foam maintains its shallow slope with 
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FIG. 8: Density profiles for piston speeds from 1 to 20 km/s, 
respectively. The snapshot times are chosen so as to capture 
the same piston location near breakout. 
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FIG. 9: Smoothed density profiles at increasing time t (in 
arbitrary units) for f/p = 5 km/s simulation. The profiles 
are centered with respect to the compaction front position 
and smoothed to demonstrate the form change. The dotted 
vertical lines show the position of the piston. 



the contact density at the piston position not changing. 
This is shown in Fig. [9] 

Mass and momentum conservation applied via the sim- 
ple snowplow model can explain the foam densification of 
the compaction front 13|. Conservation of mass is the 
essential ingredient for predicting the front speed and 
yields the simple relation 



Uf 



(1) 



where po and pc are the non-compacted foam density 
and the compacted foam density. Equation ([T]) has a 
weak dependence on the po and pc and, as such, does 
not place a high demand on modeling accuracy. Thus 



5 



40 



e 10 

Q 
LU 
LU 
Q. 

</) 



O 



experiment 
9 simulation 
(~) snowplow 



100 



10 

PISTON SPEED (km/s) 



30 



FIG. 10: Comparison of the front speeds (7/ as a function of 
the piston speeds Up for the simulations. The diamonds show 
results of experiments [l3]. The circles show the results of 
MD simulations. The hexagons show estimates based on the 
Eq. (HI). 



one should not be surprised at the agreement between 
our simple model (that assumes billiard ball atoms that 
remain neutral irrespective of temperature and pressure) 



and the very recent (but limited) experiment data 
shown in Fig. [TO] 

Also evident in Fig. [10] is good agreement between the 
simulated speeds and the snowplow model; i.e., the jump 
condition prediction of Eq. ([T|) . The compacted regions 
do not have a strictly flat (constant) density, but assum- 
ing an average density is sufficient to obtain excellent 
agreement. While this agreement is not a sensitive mea- 
sure of our knowledge of the compaction dynamics, this 
insensivity plays an important role in estimating ioniza- 
tion effects discussed later in the paper. 

Similarly, the atomic virial pressure of Kirkwood and 
the snowplow estimate. 



P = poUfUp 



(2) 



are in agreement. See Fig. [TTJ This is a consequence of 
momentum conservation across the front, a global con- 
servation requirement that is necessarily satisfied. The 
snowplow pressure P is calculated using the measured 
front speed from simulation. Again, we realize that such 
agreement is not a sensitive measure of our knowledge of 
the compaction EOS by simply examining Eq. 

We mention an important point in calculating the pres- 
sure from the virial equation. Care must be taken to use 
the local density in this expression for calculating the 
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FIG. 11: Comparison of pressures as a function of piston 
speeds for the simulations and the snowplow jump predic- 
tions. The error band is estimated using the simulation un- 
certainly in the virial pressure. 



pressure profile. Using the local density, one gets the 
correct virial pressure which is constant throughout the 
compacted foam (even though the density is not) and 
equal to the pressure estimated from the conservation of 
momentum requirement. 

The equation of state EOS of the compacted copper 
foam is presented in Fig ll2l The figure shows both non- 
ionizing and ionizing states. The applicability of the non- 
ionizing EOS for real copper foam is in doubt much be- 
yond Up = 10 km/s because of ionization of the copper 
atoms. To estimate the effect of ionization on the EOS, 
we consider two different methods. The first method uses 
tabulated EOS for hot-dense matter based on sophisti- 
cated model calculations and experiment [l5, 16 1. 

In our earlier discussion of the snowplow model, we 



emphasized that global conservation of mass and mo- 
mentum leading to Equations ^ and ([2]) demonstrate 
a weak density dependence of front speed, as well as the 
pressure. Furthermore, the front speed found from simu- 
lation agrees with the experimental measurements. This 
leads us to assume that the density and pressure for a 
given simulation is the same for both the neutral and 
ionized state, and we employed these tables to determine 
the temperatures of the ionized state. The results are 
shown in Fig. [121 We note that the ionized state temper- 
ature at Up — 20 km/s is a factor of four lower than the 
MD resuh. 

For the second method, we invoke conservation of en- 
ergy in the ideal gas (IG) approximation [13,]. We discuss 
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FIG. 12: Equation of state of compacted copper foam. The 
triangles show results from the MD simulations. Downward 
pointing triangles show results using the tables (see text for 
details), hexagons show results from an ideal gas approxima- 
tion (see text for details) . The curves are the results of fits to 
the points and are included to guide the eye. 



gas approximation, the energy is 



Eo = -NTo 
for the neutral system and 

= In{1 + Z,) T, + Ne, 



(3) 



(4) 



for the ionized system where average energy used to 
free Zi electrons from an atom. Assuming the energy 
and the number of atoms are conserved, we find then the 
following relation: 



To = (1 + Zi) + -e, 



(5) 



The ionization energies of copper are to +1: 745.5 
kJ/mol, +1 to +2: 1958 kJ/mol, +2 to +3: 3554 kJ/mol, 
+3 to +4: 5326 kJ/mol. Using the Eq (P between the 
two temperatures, we go up the ionization scale until 
it predict the lowest positive ionized temperature. In 
Fig. [TH the ideal gas (IG) predictions are plotted as 
hexagons. At the highest pressure, the IG ionization low- 
ers the temperature by an order-of-magnitude. This is a 
very simplistic model is pedagogical and shows explicitly 
why there would be a lower temperature due to ioniza- 
tion. 
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FIG. 13: This studys pressure-density EOS compared to the 
known copper Hugoniot EOS. The solid line is drawn to guide 
the eye. 
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this IG approach for pedagogical reasons even though the 
table lookup is more accurate. We denote as state 0, the 
crushed foam where the copper atoms are constrained to 
be neutral. This is our MD case. Energy, number of 
atoms and average number of electrons ionized per atom 
are denoted hy E, N and Z , respectively. In the ideal 



FIG. 14: The speed of the compaction front as a function of 
the speed of the piston. Diamonds show experimental results 
[T3 |. Downward pointing triangles show results of the MD 
simulation; the solid line is a fit to the MD results. The 
dotted line shows the theoretical copper Hugoniot [H, [l^ . 

Figure [T3] shows the pressure-density EOS from the 
MD simulations compared to the theoretical copper 
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Hugoniot [15|, |l6|. We note that the compacted foam 
state is not accurately represented by the shock hugoniot 
of low density copper foam. This is not surprising from 
all that we have learned earlier. One must be cautious 
in attempting to characterize the high-energy crushing 
of foam as a shock phenomenon. However, this does not 
preclude the very good agreement of Uf(Up) between the 
MD simulations and the copper foam hugoniot. Compar- 
ison oi Uf{Up) for MD, experiment and the copper foam 
hugoniot is shown in Fig. [Ml The agreement between all 
three is good, though MD and experiment is best. 



SUMMARY 

This study was driven by the NIF Energy and Science 
programs. We investigated the generic features of the 
compaction dynamics of metallic nano-foams as related 
to current ICF and science interests. A successful atom- 
istic method was implemented for creating a computer 
nano-foam. The microscopic processes in the compaction 
dynamics showed that a strong compaction front passing 
through a foam induces acceleration, crushing, ablation 
and mixing of the foam material within the porous back- 
ground. Depending on the piston speed, the ejecta could 
be a solid front, a liquid, an imperfect vapor, or a hot 
gas. The MD simulations agreed with limited experi- 
mental measurements. As should be expected, the sim- 
ple snowplow model gave an overall accurate macroscopic 
description of the dynamics. The EOS of the compacted 
foam is a direct consequence of the atomistic simulation, 
though its accuracy depends of a robust description of 
the atomic behavior over a wide range temperatures and 
pressures. Our simple interatomic potential did not fulfill 
that requirement because of the constrained neutrality of 
the atoms. However, simple model estimates of the ion- 
ization were made and are believed to be reasonable. The 
crushing of the foam is a highly complex evolving state 
as one moves from the front toward the piston. 
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